NMR relaxation of quantum spin chains in magnetic fields 
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We investigate NMR relaxation rates 1/Ti of quantum spin chains in magnetic fields. Universal 
properties for the divergence behavior of 1/Ti are obtained in the Tomonaga-Luttinger- liquid state. 
\£) , The results are discussed in comparison with experimental results. 
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I. INTRODUCTION 



One-dimensional (ID) quantum spin systems with an energy gap above a singlet ground state have attracted a 
great amount of attention both theoretically and experimentally. When the magnetic field is applied, the energy 
^ \ gap vanishes at the lower critical field H c \ due to the Zeeman effect and a quantum phase transition from a gapped 
spin liquid to the Tomonaga-Luttinger liquid (TLL) takes place [H 0, 0, 0, H, H 0, S S For some typical 

ID gapped spin systems, it was verified theoretically that the critical exponent of the s pin correlation function in 
+J ! H cl <H <H c2 shows characteristic field dependence in each model 0, 0, 0, II H S, 0, 0, 0, tH E El, 0, El (3 with 
i/ C 2 being the upper critical field corresponding to the saturation of the magnetization. Furthermore, it was argued 
that such characteristic field dependence of the critical exponent can be detected by NMR measurements in the TLL 
regime [1, E3]- When temperature is decreased in the TLL regime, the NMR relaxation rate 1/Ti shows divergence 
behavior and its exponent is described as a function of the critical exponent of the spin correlation function. When 
the NMR measurement is performed on the nuclei located at the sites different from the electronic spins, relaxation 
occurs through a dipolar interaction between the nuclear and electron spins. In this case, 1/Ti is expressed as a sum 
of the longitudinal and transverse relaxation processes in magnetic fields. The divergence of 1/Ti with decreasing 
I ' temperature is caused by one of these relaxation processes [ItJ ■ 
, Stimulated by these theoretical studies, NMR relaxation rates in H c i < H < H C 2 were measured in the Haldane- 
gap compound (CHQ^NNifNC^ls [lSl ] . the 5=1/2 bond-alternating spin-chain compound pentafluorophenyl nitronyl 
nitroxide (F 5 PNN) [H], and the S = 1 bond-alternating spin-chain compound Ni(CgH24N4)(N0 2 )ClC>4 [if. In these 
experiments, the increase of 1/Ti in the gapless regime was observed with decreasing temperature. However, the 
field dependence of the divergence exponent derived from the experiments is still controversial. To develop precise 
evaluation, it is indispensable to clarify the field and temperature dependences of factors other than the divergence 
' ones in the expression of 1 jT\ . 

In this paper, we investigate the NMR relaxation rate of quantum spin chains in the TLL state. In Sec. II, we 
evaluate the field and temperature dependences of factors other than the divergence ones of 1/Ti in the TLL state. 
A criterion that the temperature dependence appears only in the power-law behavior is obtained. In Sec. Ill, we 
discuss the divergence behavior of 1/Ti in connection with theoretical results obtained so far for some ID gapped 
spin systems in magnetic fields. We further discuss the experimental results in comparison with the present results. 
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II. FIELD AND TEMPERATURE DEPENDENCES OF 1/Ti IN TLL STATE 



Let us start our discussion from the formula of 1 /T\ : 

i = / p. hm l ^ [^ 2 (g)Im X ^(g,c) + F x ((,)Imx£(g, W )] , (1) 

where F zz (q) and F 1 - (q) are the longitudinal and transverse components of the hyperfine form factor, and Xr (l-> w ) 
and are the longitudinal and transverse dynamical susceptibilities defined as 

X z R z (q,u) = [dx [ At{S z {x,r)S z {Q,Q))\r^t+^ l( ^ qx \ (2) 
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xiM = \ JdxJ dt[{S + (x,T)S-(0,0)) + (S-(x,r)S + (0,O))] T ^ it+o+ ^ t -' lx) . (3) 

Spin correlation functions in the one-component TLL are expressed as 

(S z (x,t)S z (0, 0)} = m 2 + d [(a + ivr)- 2 + (x - ivr)- 2 ] + C 2 cos(2k F x)(x 2 + v 2 t 2 )- 29 • • • , (4) 



(S + (x,t)S~ (0,0)} + (S~(x,t)S + (0,0)) = C* 3 cos(7ra;)(x 2 + vV 
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(5) 



where m is the magnetization satisfying < m < 1/2, 2fcp = 7r(l — 2m), and i> and r are the velocity of the 
spinon excitation and the imaginary time, respectively. C's are positive constants, which have solely field dependence 
characteristic of the model. 

The expressions for Xr(<1^) and Xr(9? w ) a t finite temperatures were obtained in Ref. [4| by conformally mapping 
vr±ix onto the Matsubara strip as vr±ix — > (hv/irk B T) sin[(wr ± ix)irk B T /hv] in ^ and ([5]), and then performing 
Fourier transform with respect to x and t in ([2]) and ([3]). We assume that the q dependence of the hyperfine form 
factor is weaker than that of the dynamical susceptibility around the gapless point and that the hyperfine form factor 
takes the value at the wave number of the gapless point. Using the dynamical susceptibilities and this assumption, 
we obtain the following expression for 1/Ti, 
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where 27 = 2g + l/2g — 2, 5ik(q) = \q — k\ denotes the linear dispersion relation of the spinon excitation around the 
gapless point q = k (k = 2fcp,7r, and ir — 2k F ), Co — hoo / Aitk B T and q = hvq/Airk B T are the dimensionless variables, 
and 



I 7TO; / OL OL\ / OL OL\ 

n(w, Sh(q)- or) = - sin(— )B (-i [Co + 5~e k (q)} + -, 1 - - J B (-i [Q - Si k (q)} + -,!--] 



(7) 



A((D,^ 7r _ 2feF (g);2 7 ) 

= isin(7r7)|s (- i{Q + 5e n - 2 k F (?)] + ^ + 1, -1 - 7) B - ^-2k P {q)] + 2' 1_7 ) 

+ fl(-j[fl + $e w _a fcF (-g)] + ^ + 1,-1-7) s(-i[w-5e 7r _ 2fep (-g)] + |,l-7) }, 



(8) 



with -B(x, y) being the beta function. Note that the exponents of the power-law behavior about T were obtained in 
Ref. [|. 

In Fig. 1, we show the spectral weights of lmU(Co, 8ik(q); a) and ImA(d), <5e 7r _2/c F (<?); 27) for several values of a, 
where 27 = a + l/a — 2. Around the gapless point Sq — \q~ k\ ~ 0, the peak appears at Co = 0, indicating the presence 
of the overdamped spinon excitation. The NMR relaxation is dominated by the contribution of these excitations [2l| . 
As a is decreased in ImII(<I>, 5ek(q)', ot) and a approaches unity in ImA(d), 5e 7r -2k F (q)\ 2 7 ), the overdamped peak grows 
and its width becomes narrow. Note that for g = 1/2 (a = 1), where the divergence of 1/Ti vanishes, H(Co, Se^ifi); 1) 
takes the same form as the universal dynamical staggered susceptibility of an SU(2) spin chain obtained in Ref. [21j | : 



tt r(l/4-i[D+g£ x (g)])r(l/4-i[cD-g£^(g)]) 



The spectral profile of its imaginary part adequately scaled 



S) r(3/4-i[S>+i5e w (g)])r(3/4-i[£D-5e w (8)]) 1 

agrees well with each other. 

We next investigate the integrated values with respect to q in the expression ©. For the appearance of 
the TLL state, the small q region of the linear dispersion curve around the gapless point k(— 2fcp,7r, and 
7r — 2k F ) has to be considered. We thus introduce the upper cutoff of the q integral around the gapless point as 
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FIG. 2: The integrated values of It(a;T) and J„-_2fc F (27; T) as a function of 5g m ax,fc 
7r — 2&f and 2j — a + l/a — 2. 



k\, where k = 2£;F,7r, and 



h{a;T) = lim^o /j™"'* dq!mIl(Lu, Se k (q); a) and J^fep (27; T) = lim^o f-T^C ^ d( ? ImA (^> ^-2k F (q); 27)- As 
temperature is decreased, g ma x,/c becomes large. In Fig. 2, Ik(a;T) and J„-_2fc F (27; T) are shown as a function of 
fe|(= <5?fc(<7max,fc))- For small <5q maXifc , I k (a;T) and J 7r _ 2 fc F (27; T) depend on &7 max ,fc- The results 
indicate that the temperature dependence appears not only in the power-law factors but also in their coefficients. For 
<5?max,fc > 0.4, on the other hand, i/-(a;T) and J w _2fc F (27;T) take constants, which depend only on a. Since a(= 2g 
or 1/2^) has solely field dependence, I 2 k F (2g;T), I 7r (l/2g;T), and J n -2k F (27; T) are independent of temperature. In 
this case, therefore, i/ie temperature dependence of 1/Ti appears solely in the power-law behavior as 
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In Fig. 3, we show the g dependences of l2k F (2g),I- K (l/2g), and J 7r _2fc F (27). The corresponding factor 1/47T of the 
first term in §§§ is negligibly small in this scale. We find that for g > 1/2 the factor (2-nk B T jhv) 1 ^ 291 in the third 
term in ^ shows divergence and its factor I n (l/2g) becomes the largest, while for g < 1/2 the factor (2nk B T '/hv) 29 ^ 1 
in the second term shows divergence and its factor l2k F (2g) becomes the largest. The results are universal and hold 
irrespective of the model. Features of the model emerge in the field dependences of g, v and C"s in the expression ([5]). 



If we assume that the hyperfine interaction is isotropic and F zz (0) ~ F zz (2k-p) ~ F^iir) ~ F^(ir — 2ky) as usual, 



4 



125 



100 

75 
50 
25 



FIG. 3: The g dependence of I 2 fc P (2g), I n (l/2g), and J^. 2 k F {2-y), where 27 = 2g + l/2g - 2. 

the factors parenthesized by [• • •] of the four terms in ([9]) have almost the same values in a given magnetic field. Since 
the temperature dependence of l/Ti is measured in a fixed magnetic field, the factors other than those concerning 
the power-law of T in © can be regarded as constants. 

III. DISCUSSION 

We now discuss the divergence behavior of 1/Ti in connection with theoretical results obtained so far for some 
ID gapped spin systems in magnetic fields. By making use of field theoretical and numerical techniques, the field 
dependence of g was successfully obtained in several models. It was verified that g > 1/2 is satisfied in the TLL 
re gim e of the S = 1 isotropic [3, [3 fl3| and anisotropic spin chains, the S — 1/2 bond-alternating spin chain 
[313, the 5=1/2 two-leg' spin ladderjy, d, [^, [3 , an d the 5=1 bond-alternating spin chain [3- In these models, 
the transverse staggered spin correlation is dominant in H c i < H < H C 2- Therefore, the divergence behavior of the 
NMR relaxation rate is fitted well with 1/Ti = A^T 1 / 2 ^ 1 . 

On the other hand, in the 5=1/2 bond-alternating spin chain with a next-nearest- neig hbor interaction, the 
region where g < 1/2 is satisfied emerges around the half field between H c i and H c2 [H, [3 3- In other fields : 
g > 1/2 is satisfied. Such a feature is caused by the change in the dominant spin correlation in magnetic fields: 
Around H ~ (H c i + H c2 )/2 the longitudinal 2fcp spin correlation is dominant, while in other fields the transverse 
staggered spin correlation is dominant [3- In this system, accordingly, the divergence behavior of 1/T\ around 
H ~ (H cl + H c2 )/2 is fitted well with 1/Ti = A zz T 2g ~ 1 or 1/Ti = A ± T 1 / 2g ~ 1 + A ZZ T 29 ~ X . In other region of 
magnetic fields, in particular H ~ H c \ and H c2 , 1/T\ — A ± T 1 / 2g ~ 1 is adequate. 

As mentioned above, a criterion for the appearance of the temperature dependence only in the power-law divergence 
of 1/T\ is evaluated to be ^g max .fc > 0.4. We compare the temperature region derived from this criterion with that 
observed in the experiments. In the Haldane-gap compound (CH3)4NNi(NC>2)3, the divergence of 1/Ti was observed 
in 0.6K < T < 2K [3- According to the numerical calculation, the dispersion curve of the 5 = 1 Heisenberg 
model in a small magnetization region is described well as e(q) ~ 2J| sing|. The coupling constant is evaluated to 
be J = 12K [3- We estimate approximately that the linear dispersion curve may hold in the low-energy region as 
< e < 0.2 x 2 J — 4.8K. Applying 4.8K into the criterion 0.4 < 5<? maX)7r (= 5e w (g maX)W )) = 4.8K/47rT, we evaluate the 
temperature region for the divergence of 1/Ti as T < 0.95K. This temperature region overlaps with that observed in 
the experiments. Therefore, the field dependence of the divergence exponent can be well analyzed in this compound. 

In F5PNN, the divergence of 1/Ti was observed in 0.2K < T < IK 19]. From the numerical calculation, the 
dispersion curve of the corresponding model is described as e(q) ~ 0.9J|sinq| in a small magnetization region. The 
coupling constant is evaluated as J = 5.6K [23| . In the same way, the temperature region for the linear dispersion 
curve is approximately estimated as < e < 0.2 x 0.9 J ~ IK, leading to the power-law divergence region as T < 0.2K. 
This temperature region lies below that observed in the experiments. To develop more precise evaluation of the 
divergence exponent, measurements of 1/Ti in lower temperatures are necessary. 

In summary, we have investigated the NMR relaxation rate of quantum spin chains in magnetic fields. The field 
and temperature dependences of the power-law divergence of 1/Ti have been evaluated in the TLL regime. On the 
basis of the results, experimental results for some typical gapped spin chains in magnetic fields have been discussed. 
We hope that the present analyses are useful to investigate the TLL nature of quantum spin chains in magnetic fields 
via the field dependence of the power-law divergence of 1/Tl. 
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